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Abstract 

The fact that the ground-state first-order density matrix for Bosons in a parabolic 
potential with interparticle harmonic interactions is known in exact form is here 
exploited to study collective excitations in the weak-coupling regime. Oscillations 
about the ground-state density are treated analytically by a linearized equation of 
motion which includes a kinetic energy contribution. We show that the dipole mode 
has the frequency of the bare trap, in accord with the Kohn theorem, and derive 
explicit expressions for the frequencies of the higher-multipole modes in terms of a 
frequency renormalized by the interactions. 

PACS numbers: 03.75.Fi, 05.30.Jp 



1 Introduction 



In a very recent investigation, Amoruso et al. [1] have given a comparative 
discussion of collective excitations in dilute Fermi and Bose gases subject to 
harmonic confinement at zero temperature. For this purpose these workers 
employed a linearized form of the equation of motion for the density profile 
n(R, t) in the Hartree approximation [2-4], which is 
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where V ex t(R) = muo\ R 2 /2 is the confining potential, v(R — x) is the inter- 
particle interaction and Il Q , i/ g(R, t) is the kinetic stress tensor. 

Whereas for a Bose-condensed cloud Amoruso et at. [1] chose the interactions 
to be of contact form v (R — x) = g5(R — x) and adopted the Thomas- Fermi 
approximation (corresponding to n a)i g(R, t) — 0), the focus here is the N-Boson 
model discussed by earlier workers [5, 6]. In this model, all Boson properties 
are governed by harmonic forces, both confinement (which is usual) and, in 
contrast to the contact interactions in ref. [1], 

t;(R-x) = ±i 7 2 (R-x) 2 . (2) 

Notice that the positive (negative) sign in Eq. (2) corresponds to attractions 
(repulsions). 

One merit of this model is that the equilibrium density profile n (R) has the 
exact form 

n (R) = N(K N /Tr) 3/2 exp(-K N R 2 ) (3) 

with 

KN - (N-l)u ho + u N (4) 

and where 

u>% = ul ± iV 7 2 (5) 

N being the number of Bosons in the cloud. In the repulsive case one is assum- 
ing that the interactions are not so strong as to overcome the confinement. 

In this Letter we show that within this model it is possible to determine in a 
self-consistent manner the equilibrium density profile and the collective exci- 
tation spectrum of the full Eq. (1) in the weak coupling limit. The advantage 
of having analytic results is that they can be used to study and test different 
ideas and numerical methods. 



2 Equilibrium density profile at weak coupling 



We have already pointed out that Eq. (1) involves the Hartree approximation, 
corresponding to the neglect of quantal fluctuations in treating the potential 
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energy terms. We introduce a parallel treatment of the kinetic stress ten- 
sor, by using a decoupling approximation on the one-body density matrix in 
terms of the product of two condensate wave functions. Namely, we write 
(V' t (x / ,t)-?/'(x",t)) ~ $*(x',t)$(x",t) and correspondingly find 



]T V Q V^IT^R, t) = — V R • I n(R, t) V R 



VVn(R,t) 



.(6) 



In the same approximation the Boson density is given by n(R, t) = |$(R, t)\ 2 , 
i.e. we neglect the depletion of the condensate due to the interactions. This 
approximate scheme has been shown to be equivalent to adopting the Gross- 
Pitaevskii equation for the condensate [4]. 

In this approximation the equilibrium density profile n (R) is the solution of 
the following equation, 
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V ext (R) ± \ f d 3 x (R - x) 2 n (x) ^=vV^o(R) = » , (7) 

2 J 2myn (R) 



\x being the chemical potential. We choose as an Ansatz for the solution of 
Eq. (7) a gaussian profile normalized to the total number of particles: 

n (R) = N(K/7r) 3/2 exp(-KR 2 ) . (8) 



Upon substitution of Eq. (8) in Eq. (7) and equating to zero the coefficients 
of the R° and R 2 terms we find 



and 

k 2 = co 2 ho ± N^ 2 . (10) 

Of course, the width of the equilibrium profile is narrowed (broadened) by 
attractive (repulsive) interactions. 

We can now compare the approximate result in Eqs. (8) and (10) with the 
exact result in Eqs. (3) - (5). It is seen at once that the two results agree in the 
weak coupling limit (u> 2 ^> N'-f 2 ), where Eqs. (4) and (5) yield k,% ~ uj 2 ±N'~f 2 
(we are neglecting unity relative to N). This comparison emphasizes the limits 
of validity of our approach, which involves the neglect of quantum fluctuations 
and of the depletion of the condensate. 
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3 Collective excitations at weak coupling 



We proceed to evaluate the small-amplitude oscillations of the Bose cloud 
around the approximate density profile given in Eq. (8). We set n(R,t) = 
no(R) + 7ii(R, t) in the equation of motion (1) with the approximate form of 
the kinetic stress tensor in Eq. (6) and linearize it in the fluctuation ni(R, t). 
We emphasize that this procedure is treating the equilibrium profile and the 
dynamic fluctuations of the profile in a consistent manner. The linearized 
equation of motion is 



= -V R • n R V R 

at 1 m 
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Carrying out a Fourier transform with respect to the time variable and using 
Eq. (8), the equation obeyed by ni(H,u>) becomes 



cj 2 m(R, u) = ^27 2 /tn (R)R • d(u) 

+ i{v R A(R,c) + |A[ i? 3 A(R ^ ) 



;i2) 



where 



d 

A(R, uj) = VrWi(R, u) + 2KR—n 1 (R, u) + 6/cni(R, u) (13) 

oR 



and we have defined the dipole moment of the excitation as 

d(u) = J d 3 xx.n 1 (-x.,uj) . 



(14) 



Furthermore, we have used rescaled units such that m = 1 and h — 1. 

We now focus first of all on the equation of motion for a dipolar density 
fluctuation, which is obtained from Eq. (12) by taking its dipole moment. 
Setting 

ni(R,w) oc RY lm {6,<i))exp{-KR 2 ) (15) 
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in this case, it is easily seen that A(R,u) oc ni(R, uS). The average of the 
Hartree term is 

J d 3 RR a Rpn (R) = (N/2n)5 af3 (16) 



while the average of the contribution from the kinetic stress tensor yields 



2k d 



d 6 R R a I V^A(R, u) + — — [R s A(R,u)\ = k d a {uS) . (17) 



The equation of motion of the dipole therefore takes the form 

u 2 d a {u) = {K 2 TN 1 2 )d a {u) . (18) 

From Eqs. (18) and (10) we see that the dipole mode oscillates at the bare 
trap frequency, i.e. to — Uho, in agreement with the Kohn theorem [7-9]. The 
theorem would be violated if we used the exact equilibrium profile (3) - (5) in 
the approximate equation of motion (11). 

We next turn to consider the dispersion relation for all the other modes. What 
is shown below is that this dispersion behaviour for the higher-mult ipole modes 
is given by multiples of the renormalized harmonic oscillator frequency k. 

We begin by taking as an Ansatz for the density fluctuations about the spher- 
ical equilibrium profile the form 

n 1 (R,u) = $ (R)$ nlm (R) (19) 



where $ (R) oc exp(-nR 2 /2) and $ n , m (R) oc i^ +1 / 2 (/ti? 2 )r im (fl, 0) exp(-/ti? 2 /2) 
are the ground-state and excited-state wave functions of a non-interacting gas 
trapped by an harmonic (3D spherical) oscillator of frequency k, L 1 + 1 I 2 (kR 2 ) 
being the Laguerre polynomials. Consistent with this Ansatz, the Hartree term 
vanishes for all higher multipoles, from orthogonality of $ n ^ m (R) to the wave 
function R a § (R) of the dipolar excitation. We are therefore left with the 
following equation of motion for the density fluctuation: 



cj 2 m(R, u) = -V R • { n (R)V R 



=V 



2 m(R, 



vSS vWr) 
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(20) 



Since the equilibrium density is already of the form appropriate to the ground 
state of an harmonic oscillator, Eq. (20) is identically satisfied by the Ansatz 
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made in Eq. (19). Therefore, the frequencies of the modes with n > 2 are 
uj = nK, as anticipated. 



4 Summary 



In summary, we have shown that for the T = limit of the N-Boson model 
in a parabolic potential and with harmonic interparticle interactions, the lin- 
earized equation of motion involves kinetic contributions of the kind exhibited 
in Eq. (11). The corresponding Gross-Pitaevskii equation (7) for the equi- 
librium density profile has a solution which agrees with the weak-coupling 
limit of the exact profile given in ref. [5]. Inputting this equilibrium profile 
in Eq. (11) enables the dispersion relation of the collective excitations to be 
obtained analyticaly. The frequency of the dipole mode is, in accord with the 
Kohn theorem [7 - 9], that of the bare trap. The frequencies of the higher- 
multipole modes are directly related to the trap frequency renormalized by 
the interactions. 
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